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ABSTRACT 

We use three-dimensional MHD simulations with anisotropic thermal conduction to 
study turbulence due to the magnetothermal instability (MTI) in the intracluster 
medium (ICM) of galaxy clusters. The MTI grows on timescales of < 1 Gyr and is 
capable of driving vigorous, sustained turbulence in the outer parts of galaxy clusters if 
the temperature gradient is maintained in spite of the rapid thermal conduction. If this 
is the case, turbulence due to the MTI can provide up to 5-30% of the pressure support 
beyond rsoo in galaxy clusters, an effect that is strongest for hot, massive clusters. The 
turbulence driven by the MTI is generally additive to other sources of turbulence in the 
ICM, such as that produced by structure formation. This new source of non-thermal 
pressure support reduces the observed Sunyaev-Zel'dovich (SZ) signal and X-ray pres- 
sure gradient for a given cluster mass and introduces a cluster mass and temperature 
gradient-dependent bias in SZ and X-ray mass estimates of clusters. This additional 
physics may also need to be taken into account when estimating the matter power spec- 
trum normalization, ogi through simulation templates from the observed amplitude of 
the SZ power spectrum. 

Key words: convection — galaxies: clusters: intracluster medium — instabilities — 
turbulence — X-rays: galaxies: clusters 



1 INTRODUCTION 

Clusters of galaxies are on the exponential tail of the mass 
function of gravitationally bound objects in the universe. As 
a result, fundamental cosmological parameters, e.g., as, are 
very sensitive to the statistics and evolution of the cluster 
population. Unfortunately, only ~ 12% of the total mass of 
a cluster is contained in the hot intracluster medium (ICM) 
that is readily observable. Thus for many observational probes 
the total cluster mass must be estimated by either fitting 
to hydrostatic equilibrium or by assuming scaling relations 
between total mass M and other observables, e.g., M-L x , 
where L x is the total X-ray luminosity. 

A new window for studying galaxy clusters has recently 
been opened with measurements of the Sunyaev-Zel'dovich 
(SZ) effect, which is produced by the inverse Compton scat- 
tering of CMB photons off the hot ICM. Recently a vari- 
ety of surveys, including the South Pole Telescope (SPT) 
and the Atacama Cosmology Telescope (A CT) have compiled 
the fi rst SZ-selected cluster samples (e.g., IVanderlinde et all 
120101 ) ■ The SZ signal scales linearly with gas density (pro- 
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portional to total pressure) allowing observations to explore 
larger cluster radii than is typically possible with X-ray obser- 
vations. A disadvantage of moving to larger radii is that one is 
potentially more susceptible to the turbulence and bulk flows 
that result from structure formation processes. These sources 
of non-thermal pressure support may significantly complicate 
the extraction of cosmological parameters from SZ surveys 
(|Shaw et al.ll20ld) . 

We have learned in recent years that convection in the 
ICM is very different from the more familiar convection in 
stars. In a dilute, magnetized medium like the ICM, the 
mean free path is much larger than the gyroradius and there- 
fore thermal conduction is entirely anisotropic along mag- 
netic fi eld lines. In this regime the magnetot hermal instability 
(MTI; |BaIbuj|2OO0l : IParrish fc Stone! 120051 ) can tap into the 
radially decreasing temperature gradient to drive convection 
regardless of the background entropy gradient. The MTI in 
Cartesian simulations with a fixed temperature gradient is ca- 
pable of driving vigorous convection (approaching supersonic 
velocities) and a mag netic dynamo that is similar to that of 
adiabatic convection (|McCourt et al.ll201ll ). 

IParrish et"afl (120081 ) previously studied the MTI in global 
galaxy cluster models. In their simulations, the rapid thermal 
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conduction made the outer parts of the cluster isothermal af- 
ter ~ 1 Gyr, suppressing the free energy driving the MTI. 
However, there are observational studies that suggest that 
most clusters have temperature profiles that dec line with ra- 
dius outside a few hundred kpc jPratt et alj|2007h . Motivated 
by these results, we revisit 3D MHD simulations of the MTI 
in the ICM; in contrast to previous work, we now prescribe 
a fixed temperature gradient in the outskirts of clusters. In 
analyzing the results of these simulations, we focus on quan- 
tifying the non-thermal pressure support that results from 
MTI-driven turbulence. In Section 2, we describe our method- 
ology and fiducial cluster model. In Section 3, we describe the 
results of our numerical experiments; we discuss their impli- 
cations for SZ and X-ray studies of galaxy clusters in Section 
4. 



2 METHOD AND MODELS 

We solve the usual equations of magnetohydrodynamics 
(MHD) with the addition of anisotropic thermal conduction. 
The MHD equations in conservative form are 
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where the symbols have their usual meaning. The total energy 
E is given by 

B B 



_ v ■ v 

E = e + p— — + 



8tt 



(5) 



where e = p/(j — 1). Throughout this paper, we assume 7 = 
5/3. The anisotropic electron heat flux is given by 

Q = -K Sp bb ■ VT, (6) 

where ks p is the Spitzer conductivity (|Spitzerlll962h and b is 
a unit vector in the direction of the magnetic field. 

We simulate the fully 3D time-dependent evolution of 
our model galaxy clusters using the MHD code Athena 
(|Stone et al.l [2008 ) with the addition of a module for 
anisotropic thermal conduction along magnetic field lines 
(|Parrish fc Stonell2005l : ISh arma fc Hammettll2007l V Our ini- 
tial condition is a spherically-symmetric, hot, massive clus- 
ter that roughly resembles Abell 1576 (although not fit to 
the exact parameters) with a mass of 1.6 x 10 15 M Q . We 
use a softened NFW profile with a scale radius of r s = 600 
kpc and a softening radius of 70 kpc. We initialize an at- 
mosphere in hydrostatic equilibrium using the entropy power 
law in the ACCEPT database for Abell 1576: a central en- 
tropy Ko = 186 keV cm 2 , K\ = 98 keVcm 2 , and power- law 
exponent, a = 1.38 (|Cavagnolo et all 120091 ). We use a mean 
molecular weight p ~ 0.62 which corresponds to a metallic- 
ity of approximately 1/3 solar. If we assume that our fidu- 
cial cluster is located at z = 0.1, then for the appropriate 



WMAP5 cosmological parameters rsoo = 1-09 Mpc, and the 
virial radius is r-2oo = 1.6 Mpc, where ta corresponds to an 
overdensity of A times the critical density. We do not include 
cooling, as we focus on the portion of the cluster that is well 
outside the cooling radius. 

The simulations are carried out on a (196) 3 Cartesian 
grid in a computational domain that extends from the center 
of the cluster out to ±1300 (2400) kpc for our fiducial (or 
larger domain) simulations. Within this Cartesian domain, 
we define a spherical subvolume with a radius of 1225 (2200) 
kpc from which we extract cluster properties, thus avoiding 
boundary condition effects. In this volume we initialize tan- 
gled magnetic fields with (\B\) = 10~ 8 G (pla sma ~ 10 4 - 
10 6 ) a nd a Kolmogorov power spectrum (see IParrish et alj 

(2009) for details). We have also performed runs at (96) 3 and 
(288) 3 and find that the kinetic energies are reasonably well 
converged at our fiducial resolution. The magnetic field am- 
plification is not quite as well converged, and higher resolu- 
tion simulations have higher final magnetic field strengths. 
We choose a smaller magnetic field strength than that ob- 
served in z — clusters as a guess of what the magnetic field 
strength was when clusters first formed. In order to simu- 
late clusters with a negative radial temperature gradient, we 
fix the temperature at the peak of the cluster temperature 
profile (approximately 200 kpc) and at the maximum radius 
of our model cluster to the initially-computed temperature 
values. Thus, we are imposing Dirichlet boundary conditions 
on the temperature profile. This fixed temperature gradient 
then drives the continued evolution of the MTI. Our fiducial 
model has a peak temperature of 10 keV at 200 kpcs and 
a temperature gradient of 3.5 keV (i.e., AT/T ~ 1/3) over 
~ 1 Mpc. This fixed temperature boundary co ndition repre- 
sents the key difference between this work and lParrish et all 
(|2008h . 

We also examine the interplay between other sources 
of turbulence and the MTI. Such turbulence could arise 
from galaxy wakes or structure formation, and i s almost cer- 
tainly present in the outer regions of clusters (|Nagai et alj 
l2007f > . We drive turbulence with an outer scale of 200 kpc 
and a steep spectrum (vk ~ &~ 3 ) with an energy injec- 
tion rate per unit volume such that the turbulence has a 
Mach number, M ~ 0.13, within 1 Mpc (more specifically 
e = 1.25 x 10 _30 erg cm _3 s _1 ). This turbulence is driven in 
Fourier space with random phases with a radiall y constant en- 
ergy i njection rate using the method detailed in lParrish et al.l 

(2010) . Our toy model for turbulence results in Mach num- 
bers somewhat smaller than fou nd in non-condu cting struc- 
ture formation calculations fe.g. lLau et alj[2009h . 



3 RESULTS 

Our und erstanding of th e satu ration of the MTI was eluci- 
dated in iMcCourt et al.1 f|201lh ■ The MTI is most unstable 
for dynamically weak (/3 ^> 1) magnetic fields that are per- 
pendicular to gravity. As the instability develops, the mag- 
netic field initially becomes progressively more aligned with 
the local gravitational field. The initial intuition was that the 
MTI behaved as a dual to the heat-flux-driven buoyan cy in- 
stability fHBI; lQuataertl2008l : |Parrish fc Quataertfcood ). and 
saturated by turning horizontal (_L g) magnetic field lines into 
vertical (|| g) field lines, a bulk reorientation of the magnetic 
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Figure 1. Azimuthally-averaged rms Mach number profiles for our 
fiducial cluster model (AT = 3.5 keV, r 50 o = 1.1 Mpc). The blue 
dotted line shows a simulations with anisotropic thermal conduc- 
tion and the MTI. The red dashed lines shows the same cluster 
model with no conduction (and no MTI), but with externally- 
driven turbulence. The black solid line shows simulations with both 
anisotropic conduction and externally-driven turbulence; the tur- 
bulence resulting from the MTI and the driven turbulence add 
approximately linearly in energy (quadrature in M). 



field geometry. However, if the field lines are predominantly 
vertical, horizontal (_L g) motions have essentially no restor- 
ing force (for weak magnetic fields), and any small perturba- 
tion produces a horizontal field geometry that is again MTI 
unstable^ Thus, the MTI cannot be suppressed by simply 
rearranging the magnetic field. Instead, in Cartesian simula- 
tions of approximately one scale height in size, the MTI drives 
sustained convection with rms Mach numbers of ~0.3. 

To understand the saturation of the MTI in global clus- 
ters with a fixed temperature gradient, we first consider our 
fiducial cluster model with the domain size of 1225 kpc. The 
MTI growth time is given by 

* MTI= ( 5_ dV) ~ 600 Myr. (7) 

If we assume that major mergers occur on a 5 Gyr timescale 
then there are over 8 e-folding times available for the MTI 
to grow; thus there is sufficient time for the MTI to be- 
come highly nonlinear. In approximately 2-3 Gyr, the tur- 
bulence reaches a statistical steady state. Figure [JJ shows an 
azimuthally-averaged radial profile of the Mach number in 
this steady state for our fiducial cluster model. The MTI 
drives large turbulent velocities with mean Mach numbers 
~ 0.1 — 0.2 within rsoo- We also simulated our fiducial cluster 



1 Buoyant overstabilitie s, such as those highlighted in 
iBalbus fc Reynolds! tepid) , can also be present; however, their 
effect is subdominant as they grow much more slowly than the 
MTI. 
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Figure 2. Azimuthally-averaged rms Mach number profiles for 
our fiducial cluster model with different magnitudes of the tem- 
perature gradient. The temperature gradients are labeled by the 
temperature drop from the temperature maximum at ~ 200 kpc 
to rsoo — IT Mpc. The turbulent velocities scale roughly as 
v ~ (AT) 1 / 2 . 



with anisotropic Braginskii viscosity and found that the tur- 
bulent velocities are the same to within ~5%. This negligible 
impact of viscosity on the e volution of t he MTI is consistent 
with theoretical predictions |Kunzll201ll ). 

It is also important to assess how the MTI interacts with 
other sources of turbulence — turbulence that could be driven 
by structure formation or other processes. Figure [JJ shows the 
saturated Mach number profiles for numerical experiments in 
which we add turbulent external driving to our fiducial model 
both with and without anisotropic conduction. We find that 
the externally-driven turbulence adds to the MTI-driven tur- 
bulence such that the turbulent energies roughly add linearly. 
This property corresponds to turbulent velocities and Mach 
numbers that add in quadrature: 

Mot « {M 2 MTl + M 2 tulb ) 1/2 . (8) 

We find that this scaling holds for levels of external turbu- 
lent driving both larger and smaller than shown in Figure [TJ 
provided that the Mach numbers are reasonably subsonic. 

Although the kinetic energy quickly saturates to a rela- 
tively steady state value for both the MTI-only and the MTI 
plus turbulence simulations, the growth of the magnetic field 
is comparatively slow. More quantitatively, the logarithmic 
growth rate of magnetic energy for our fiducial case with both 
conduction and driven turbulence is dln(_B 2 )/di ~ 0.4 Gyr" 1 . 
This is ~ 4 times slower than the linear growth of the MTI. 
Initially the magnetic fields are completely tangled; thus, the 
orientation is statistically isotropic: (\b ■ f|) = 0.5. For the 
MTI-only run, the magnetic geometry develops a small radial 
bias, peaking around (\b-f\) ~ 0.65 at 2 Gyr before declining 
to isotropy again. For the run with the MTI and driven turbu- 
lence, the magnetic field only briefly deviates from statistical 
isotropy, fluctuating around (|b ■ r|) = 0.5. 
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Previously iRuszkowski et all (|2010f ) performed a cos- 
mological study of the evolution of a single cluster with 
anisotropic thermal conduction. They found evidence for a 
small radial bias in the magnetic field; however, in the cos- 
mological simulations the effects of the MTI and radial infall 
were not disentangled. In light of our simulations, it is likely 
that the radial magnetic fields in the cosmological simulation 
are due to infalling substructure dragging magnetic field lines 
out. The isolated non-cosmological simulations presented here 
are necessary to separately understand these effects. Our re- 
sults also sugg est that the radial m agnetic field bias observed 
in Virgo l|Pfrommer fc Du rsi 2010) is also likely due to infall 
or galaxy motions rather than the MTI. 

For a given cluster mass, the saturated state of the MTI 
depends primarily on dlnT/dlnr, since the free energy of the 
temperature gradient drives the MTI. We adjust our initial 
equilibrium to several different temperature gradients by in- 
creasing the temperature maximum to study this effect, hold- 
ing all of the remaining model parameters fixed. Figure [2] 
shows that there is a clear trend in the resulting turbulence 
driven by the MTI, with larger temperature gradients pro- 
ducing stronger turbulence. If we calculate a power-law fit 
to the turbulence induced by the MTI at rsoo we find that 
Ai oc (AT) 0,5 , where AT is the temperature drop across 1 
Mpc. This scaling is consistent with the simple mixing length 
estimate of v 2 oc AT. 

Much of the contribution of clusters to the SZ power 
spectrum comes from scales comparable to or larger than the 
virial radius of groups and clusters. In order to understand 
the possible effects of the MTI on such SZ power spectrum 
measurements, we have extended our simulation domain to a 
larger size (~2rsoo). In these larger domain simulations, the 
MTI drives even more vigorous convection. Figure[3]shows the 
steady-state turbulent pressure profiles for our large domain 
simulations with the MTI alone, externally-driven turbulence 
alone, and both together. The MTI is able to drive convection 
that yields significant non-thermal pressure support, reach- 
ing ~ 35% of the thermal pressure near 2r 5 oo- The turbu- 
lent Mach numbers increase significantly towards the cluster 
outskirts. Our interpretation is that large convective motions 
deep in the core overshoot and continue with large momentum 
to the low-density cluster outskirts; this effe ct was also seen 
in ou r earlier stratified Cartesian simulations l|McCourt et al.l 

Hoiil). 

More quantitatively, in a shell centered at 1800 kpc of 
width 50 kpc, we find that the rms Mach number with la fluc- 
tuations is (M) ~ 0.54 ± 0.19. The maximum Mach number 
in this shell reaches 1.26. The behavior of the simulation with 
both the MTI and turbulence simultaneously is particularly 
interesting — within approximately 1200 kpc, the MTI-driven 
turbulence and external turbulence add in quadrature as seen 
previously in Figure Q] however, beyond this radius the net 
turbulence is approaching trans-sonic and the MTI is unable 
to grow effectively. If we perform the same experiment with 
the higher temperature gradient of AT = 5.5 keV, we find 
that the Mach number can reach values as high as M ~ 2.5, 
also likely due to overshoot. Unfortunately, we are not able 
to run these simulations to a full steady-state as the MHD 
integrator in Athena crashes in high Mach number, highly 
stratified, turbulence. 
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Figure 3. Azimuthally-averaged, steady-state turbulent pressure 
profiles (normalized to the thermal pressure) for the larger volume 
simulations of our fiducial cluster model (AT = 3.5 keV, r^oa = 1.1 
Mpc) . The legend is the same as Figure [l] At r > rsoo , the turbu- 
lent pressure support can be a significant fraction of the thermal 
pressure. 



4 DISCUSSION AND IMPLICATIONS 

Galaxy clusters have the potential to serve as powerful probes 
of structure formation by breaking the <j$-Qm degeneracy 
and constraining the dark energy equation of state, w. More- 
over, much of the constraining power of galaxy clusters is 
highly complementary to constraints from the CMB and other 
methods. However, for both SZ and X-ray observations, these 
applications require accurately determining the total clus- 
ter mass either directly through reconstructing a hydrostatic 
mass or indirectly through calibrated scaling relations (note 
that this is, of course, not true for lensing mass determina- 
tions). It is critical to understand any systematic errors or 
biases in these cluster mass determinations. 

Cosmological hydrodynamic simulations have shown for 
some time that turbulence produced by structure formation 
can provide significant non-thermal pressure support in clus- 
ters, particularly at large radii (e.g-. lEvra rd 1990; Vazz a et alj 
2009). Our results demonstrate that the MTI, a convective 
instability triggered by anisotropic thermal conduction along 
magnetic field lines, can also produce significant (near-sonic) 
turbulent motions in the outer parts (> 200 kpc) of clusters. 
This turbulence is the strongest at large radii ~ V2oo (Fig. [3]), 
in part because of convective overshooting from smaller radii. 
The non-thermal pressure support produced by the MTI is 
not directly accounted for in current attempts to estimate 
cluster masses from observables such as the X-ray luminosity 
or SZ signal. 

Mass measurements of galaxy clusters using Chandra or 
XMM-Newton have typically been confined to within rsoo 
where the x-ray surface brightness is reasonably large. Even 
within this radius, we find non-neglible turbulent pressure 
support ~ 3-10% (Fig. due to the MTI, with the exact 
value depending on the cluster temperature gradient. X-ray 
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observations of clusters have be en used to try to c onstrain the 
properties of dark energy (e.g.. lAllen et al.1 l2008): our results 
suggest the need for somewhat larger priors on the fraction 
of non-thermal ICM pressure support in such modeling. 

The thermal SZ effect is often characterized by the 
Compton-y parameter, which measures the fractional change 
of the CMB temperature and is oc n e T e . The total SZ signal 
is then simply proportional to the total thermal pressure of 
the cluster plasma. As a result, the SZ signal correlates well 
with the total cl uster mass i n cosmological hydrodynamics 
simulations (e.g.. iNagal l2006h . Turbulence due to the MTI 
would decrease the SZ signal for a given cluster mass rela- 
tive to these scaling relations; it would also likely introduce 
additional scatter into such relations. 

In addition to finding individual clusters, SZ experi- 
ments survey the entire sky to produce an angular power 
spectrum whos e amplitude is very sensi tive to as, scaling as 
Ci oc al(Q b h) 2 l|Komatsu fc SeliaklfeOQ^ ') . The connection be- 
tween as and the SZ power spectrum requires calibration by 
simulations, and the first SZ results found as = 0.746 ±0.017, 
a v alue that was stati s tically in tension w ith other probes of 
o- 8 i|Lueker et al.ll2010l ; ISengal et al.ll201ll ). Resolving this dis- 
crepancy requires decreasing the predicted SZ power by a fac- 
tor of ~2. Both the ACT and SPT groups have attempted to 
remedy this model-specific deficiency by r unning cosmologi- 
cal simulations with additional feedback llShaw et alj 2010l; 
iBattaglia et al.l |2010| ). Analytical modeling by Ishaw et all 
<|2010h shows that increasing the fraction of non-thermal 
pressure support in galaxy clusters can significantly reduce 
SZ power at large angular scales. Thus, the additional non- 
thermal pressure support produced by the MTI at large radii 
in galaxy clusters (Fig. [3} may significantly modify the deter- 
mination of as from the SZ power spectrum. 

In Figure [2] we have shown that there is a direct con- 
nection between the level of turbulent non-thermal pressure 
support due to the MTI and the cluster temperature gradi- 
ent at large radii. This is a simple consequence of the fact 
that it is the temperature gradient that drives the instability 
in the first place. This result is particularly interesting be- 
cause it implies that constraints on the temperature gradient 
could be used to estimate the MTI contribution to the tur- 
bulent pressure, which would be very valuable in improving 
cluster mass estimates. Absent such observational constraints, 
however, Figure [5] implies that variations in cluster temper- 
ature gradients as a function of mass and/or redshift could 
introduce a subtle bias in the X-ray and SZ signals used to 
constrain cosmological parameters. In addition, in lower mass 
clusters or groups with lower virial temperatures, the thermal 
conduction time across the cluster outskirts will be longer 
than in the massive cluster model considered in this work. 
This is likely to cause the turbulent pressure support induced 
by the MTI to depend explicitly on cluster mass. Recent X- 
ray observations of Perseus with the Suzaku telescope have 
shown eviden ce that suggests clump ing of baryons near the 
virial radius l|Simionescu et al]|201ll ). Turbulence can drive 
fractional density enhancements (Sp/p) that scale as either 
M 2 for compressive fluctuations or linearly as M if the fluc- 
tuations are behaving as a passive scalar in the stratified at- 
mosphere. While not large enough in magnitude alone, the 
turbulence we find can contribute to these density enhance- 
ments. 

In our calculations we have fixed the temperature gra- 



dient at large radii in our model clusters by using Dirichlet 
boundary conditions. In reality, the temperature gradient is 
likely set by an interplay between structure formation (e.g., 
the virial shock) and thermal conduction itself. In previous 
simulations of isolated clusters with anisotropic thermal con- 
duction that did not fix the temperature gradient, thermal 
conduction was effective enough to wipe out the temperature 
gradient after ~ 1 Gyr, strong ly suppressing the le vel of turbu- 
lence generated by the MTI i|Parrish et alJl200Sft . The obser- 
vational evidence for non-zero tempera ture gradients at large 
radii in clusters (e.g.. |Pratt et al]|2007l ) motivated our choice 
of boundary conditions in this work. Physically, we regard this 
choice as a proxy for cosmological physics not included in our 
simulations (e.g., heating by substructure moving through the 
ICM and/or compressional heating and inflow after the virial 
shock). It remains to be seen, however, whether significant 
temperature gradients can indeed be maintained. The only 
cosmological simulation with anisotropic thermal conduction 
to date found dT/dr < at large radii, with the magnitude 
of the temperature gra dient comparable to th at assumed in 
this work (see Fig. 2 of [Ruszkows ki et al.ll2Q10h . 

The non-cosmological simulations presented in this paper 
are not sufficient to accurately determine the contribution of 
the MTI to the turbulent pressure support in clusters with- 
out calibration by comparison to cosmological simulations. 
This is particularly true at large radii ~ r2oo where our iso- 
lated cluster simulations predict near sonic motions for even 
modest temperature gradients (Fig. [3j . We suspect that the 
convective overshooting that produces these large velocities 
at large radii is likely to be generic. In the future, a suite of 
cosmological simulations with anisotropic thermal conduction 
will be necessary to fully assess the effects of the MTI on SZ 
and X-ray observations of clusters. Such simulations are now 
possible but isolated cluster simulations like those presented 
here are still necessary to pinpoint and understand the key 
physics. 
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